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ABSTRACT 

Whether analytic exact vacuum(electrovacuum) solutions of the Einstein(Einstein-Maxwell) field equations 
can accurately describe or not the exterior spacetime of compact stars remains still an interesting open question 
in Relativistic Astrophysics. As an attempt to establish their level of accuracy, the radii of the Innermost Stable 
Circular Orbits (ISCOs) of test particles given by analytic exterior spacetime geometries have been compared 
with the ones given by numerical solutions for neutron stars (NSs) obeying a realistic equation of state (EoS). 
It has been so shown that the six-parametric solution of Pachon, Rueda, and Sanabria (2006) (hereafter PRS) 
is more accurate to describe the NS ISCO radii than other analytic models. We propose here an additional 
test of accuracy for analytic exterior geometries based on the comparison of orbital frequencies of neutral test 
particles. We compute the Keplerian, frame-dragging, as well as the precession and oscillation frequencies of 
the radial and vertical motions of neutral test particles for the Kerr and PRS geometries; then we compare them 
with the numerical values obtained by Morsink and Stella (1999) for realistic NSs. We identify the role of high- 
order multipole moments such as the mass quadrupole and current octupole in the determination of the orbital 
frequencies especially in the rapid rotation regime. The results of this work are relevant to cast a separatrix 
between black hole (BH) and NS signatures as well as probe the nuclear matter EoS and NS parameters from 
the Quasi-Periodic Oscillations (QPOs) observed in Low Mass X-Ray Binaries. 
Subject headings: Exact solutions Einstein-Maxwell equations - Relativistic precession frequencies 



1. INTRODUCTION 

One of the greatest challenges of the general theory of re- 
lativity has been the construction of solutions to the Einstein- 
Maxwell field equations representing the gravitational field 
of compact stars such as neutron stars (NSs). Stationary ax- 
ially symmetric spacetimes satisfy basic properties one ex- 
pects for rotating objects, namely time symmetry and re- 
flection symmetry with respect to the rotation axis (see e.g. 
Pachon and Sanabria-Gomez 2006). The simplest stationary 
axially symmetric exact exterior vacuum solution describing 
a rotating configuration is the well-known Kerr metric (Kerr 
1963). The Kerr metric is fully described by two free pa- 
rameters: the mass M and the angular momentum J of the 
object. However, it is known from numerical models that 
the quadrupole moment of rotating NSs deviates considerably 
from the one given by the Kerr solution Qkbh = -J~/iMc~) 
(see e.g. Laarakkers and Poisson 1999, for details). 

In the mean time, a considerable number of analytic ex- 
terior solutions with a more complex multipolar structure 
than the one of the Kerr solution have been developed (see 
e.g. Manko et al. 1995, 2000; Stephani et al. 2003). Whether 
analytic exterior solutions are accurate or not to describe 
the gravitational field of compact stars is an interesting and 
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very active topic of research (see e.g. Stute and Camenzind 
2002; Berti and Stergioulas 2004; Pachon et al. 2006, and ref- 
erences therein). 

The accuracy of analytic solutions to describe the exterior 
geometry of a realistic rotating compact star has been tested 
by comparing physical properties, e.g. the radius of the Inner- 
most Stable Circular Orbit (ISCO) on the equatorial plane and 
the gravitational redshift (see Sibgatullin and Sunyaev 1998; 
Berti and Stergioulas 2004; Pachon et al. 2006, for details). 
In order to do such a comparison, the free parameters (i.e. the 
lowest multipole moments) of the analytic exterior spacetime, 
are fixed to the corresponding lowest multipole moments 
given by numerical interior solutions of the Einstein equa- 
tions, for NS realistic models (see e.g. Berti and Stergioulas 
2004). 

Following such a procedure, the solution of Manko et al. 
(2000) has been compared by Stute and Camenzind (2002) 
and by Berti and Stergioulas (2004) with the numerical 
solutions for NSs calculated by CookefaL (1994) and 
with those derived by Berti and Stergioulas (2004), respec- 
tively. However, being a generalization of the solution of 
Tomimatsu and Sato (1972), it cannot describe slowly ro- 
tating compact stars (see e.g. Berti and Stergioulas 2004), 
but the dynamics of astrophysical objects with anisotropic 
stresses (see Dubeibe et al. 2007, for details). 

Following a similar procedure, based on tests of the ISCOs 
radii on the equatorial plane of the rotating neutron stars ob- 
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tained by Berti and Stergioulas (2004), it has been shown that 
the six-parametric solution of Pachon et al. (2006) (hereafter 
PRS solution, see Sec. 2 for details) is more accurate than the 
model of Manko et al. (2000). In addition, being a generaliza- 
tion of the Kerr solution, this solution can be used for arbitrary 
rotation rates. 

Besides the ISCOs radii, there are additional physical prop- 
erties that can be computed with analytic and numerical mod- 
els and thus useful to compare and contrast the accuracy of 
analytic exact models. The aim of this article is to analyze 
the properties of orbital frequencies of neutral test particles in 
the PRS and in the Kerr geometries with especial focus on the 
Keplerian vk, frame-dragging (Lense-Thirring) vlt, as well as 
the precession(oscillation) frequencies of the radial and verti- 
cal motions, Vp(v°^) and v?(v? ), respectively. 

The relevance of these frequencies relies on the fact that 
they are often invoked to explain the Quasi-Periodic Os- 
cillations (QPOs) observed in some relativistic astrophysi- 
cal systems such as Low Mass X-ray Binaries (LMXBs), 
binary systems harboring either a NS or a black hole 
(BH) accreting matter from a companion star. For in- 
stance, within the Relativistic Precession Model (RPM) intro- 
duced by Stella and Vietri (1998); Morsink and Stella (1999); 
Stella effl/. (1999); Stella and Vietri (1999), the kHz QPOs 
are interpreted as a direct manifestation of the modes of rela- 
tivistic epicyclic motion of blobs arising at various radii r in 
the inner parts of the accretion disk around the compact object 
(see Sec. 6, for details). 

In addition to the RPM, the Keplerian, precession and os- 
cillation frequencies are used in other QPO theoretical mod- 
els (see e.g. Lin etal. 2011, for a recent comparison of the 
existing models). Due to the influence of general relativistic 
effects in the determination of such frequencies, an observa- 
tional confirmation of any of the models might lead to an out- 
standing test of general relativity in the strong field regime. 
In this line, it is of interest to compare and contrast the or- 
bital frequencies given by the Kerr solution and by the PRS 
solution (see Sec. 3), which help to establish the differences 
between possible BH and NS signatures. We emphasize in 
this article the major role of the quadrupole moment as well 
as of the octupole moment of the object, whose possible mea- 
surement can be used as a tool to test the no-hair theorem of 
BHs (see e.g. Johannsen and Psaltis 201 1) and to discriminate 
between the different theoretical models proposed to explain 
the physics at interior and exterior of the Neutron Stars. Ad- 
ditionally, in the case of NSs, the interpretation of QPOs as 
the manifestation of orbital motion frequencies might lead to 
crucial information of the NS parameters such as mass, an- 
gular momentum (see e.g. Stella and Vietri 1998; Torok et al. 
2010), and quadrupole moment (see e.g. Morsink and Stella 
1999). These parameters reveal, at the same time, invaluable 
information about the EoS of nuclear matter 

The article is organized as follows. In Sec. 2 we recall the 
properties of the PRS solution. The computation of the or- 
bital frequencies as well as the comparison of their features 
in the Kerr and in the PRS spacetimes, is shown in Sec. 3. In 
Sec. 4 we study the accuracy of the analytic formulas of the 
periastron and nodal frequencies derived by Ryan (1995) for 
stationary axially symmetric spacetimes. In Sections 5 and 
6 we discuss the accuracy of the PRS solution in describing 
the frequencies of realistic NS models and its relevance in the 
Relativistic Precession Model, respectively. The conclusions 
of this work and a discussion on possible additional effects to 



be accounted for in the determination of the orbital frequen- 
cies, e.g. the effect of magnetic dipole moment, are outlined 
in Sec. 7. 

2. THE PRS ANALYTIC EXACT SOLUTION 

We first recall the PRS analytic model (Pachon et al. 2006), 
for the exterior gravitational field of a compact object"*. In the 
stationary axisymmetric case, the simplest form of the metric 
can be written as (Papapetrou 1953) 

ds' = -f(dt - cod(pf + [e^^idp^ + dz^) + p^d4>^] , (1) 

where /, o) and y are functions of the quasi-cylindrical Weyl 
coordinates (p, z). Thus, the components of the metric tensor 
gfiv are 



f{p,z) 
git = -f(fi,z), 
gt<P = /(P> z)u}(p, z), 

p2r(p,z) 



f(p, z)oj(p, z) 



1 



gzz = gpp 



fip, Z) g^ gi 
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Using the above line element, the Einstein-Maxwell equa- 
tions can be reformulated, via Ernst's procedure in terms of 
two complex potentials &{p,z) and <l>(p, z) (Ernst 1968a,b). 
By means of Sibgatullin's integral method (Sibgatullin 1991; 
Manko and Sibgatullin 1993) this system of equations can be 
solved va 



do- e(^)iu(o-) 



vr 



/do 



dcrm^iicr) 



vr 



(6) 



(7) 



where e(z) := S(z,p = 0) and f(z) := <l'(z,p = 0). The un 
known function //(cr) must satisfy the singular integral equa 
tion 

^' fi(cr)[e(0 + ~e(rf) + imfmdo- 



£ 



(o- - r) Vl -cr2 
and the normalizing condition 



= 



X 



fi(o-)do- 



-1 vr 



(8) 



(9) 



where ^ - z + ipo; rj - z + ipr, p and z being the Weyl- 
Papapetrou quasi-cylindrical coordinates, cr, r e [-1,1], 

e(ri) := e(^), /(rj) := /(jj) and the overbar stands for com- 
plex conjugation. In (Pachon et al. 2006), the Ernst potentials 
were chosen as 



e(z) 



Z'^ - z^(m + ia) - kz + is 
z-^ + z^(m — ia) — kz + is 

qz^ + ifiz 
Z'' + z^(m — ia) — kz + is 



(10) 



* Mathematica 8.0 scripts with the solution, some limiting cases as 
well as the the calculations presented in this paper are available at 
http://www.chem.utoronto.ca/~lpachon/scripts/nstars 
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We calculate the multipole moments following the proce- 
dure of Hoenselaers and Perjes (1990). We denote the mass 
multipoles by M, while, the current (rotation) multipoles, by 
Si. The electric multipoles are denoted by Q, and the mag- 
netic ones by S,. Thus, for the PRS solution we have 

Mq - m , M2 - mk - ma^ , 
S I = ma , S3 — -ma^ + 2mak - ms , 
Qo - q, Qi- -a'q -a^ + kq, ... 



(11) 



^\ — fj. + aq , S2 — —a^n + fik — a^q + 2akq — qs , 



(12) 



This allows us to identify m as the total mass, a as the total 
angular moment per unit mass {a — J/m, being J the total 
angular moment); while k, s, q and /i are associated to the 
mass-quadrupole moment M2, current octupole 53, electric 
charge and magnetic dipole, respectively. 

The potentials (10) can be written in an alternative way, we 
mean 



3 



fi 



(13) 



with 



iu + dBj 



(14) 

, i,k + i. (15) 

iPj-mPj-Pi) 
Then, using Eqs. (6) and (10), we obtain the Ernst potentials 



£■ ^ + ^ ^ C 
6 = , O 



A-B A-B 
and the metric functions in the whole spacetime 



(16) 



^ AA-BB + CC AA-BB + CC 
f = — ^TTT — ^ ' ^ = : ' (17) 



(A - B)(A - B) 



KKY]r„ 



Im[(A + B)H - (A + B)G - CI] 
AA-BB + CC 



(18) 



where the functions A, B, C, H, G, K, and / can be found in 
the Appendix A. 

The PRS electrovacuum exact solution belongs to the ex- 
tended A^-soliton solution of the Einstein-Maxwell equations 
derived by Ruiz et al. (1995), in the particular case N -3. In 
addition, the functional form of the metric functions resem- 
bles the one derived previously by Breton et al. (1999). Be- 
sides the limiting cases discussed in Pachon et al. (2006) it is 
worth mentioning that, in the vacuum case ^ = and // = 0, 
for s - Q this solution reduces to the solution of Manko et al. 
(1995) under the same physical conditions, namely q - Q, 
c = and /7 = in Manko et al. (1995). 

3. ORBITAL MOTION FREQUENCIES ON THE EQUATORIAL PLANE 

Although for the case of compact stars contributions from 
the magnetic field could be relevant (see e.g. Bakala et al. 
2010; Sanabria-Gomez et al. 2010; Bakala et al. 2012), we fo- 
cus in this work on the frequencies of neutral particles orbit- 
ing a neutral compact object. We calculate here the Keple- 
rian vk - n^/iln), frame-dragging (Lense-Thirring) vlt - 



D.u:/{2n), radial oscillation and precession, = Q°^/(27r) 
and = Slp/{2n), and vertical oscillation and precession fre- 



quencies. 



„os 



OS 



/(2n) and = D.pl(2n), respectively. 



The geodesic motion of test particles along the radial co- 
ordinate, on the equatorial plane z = 0, is governed by the 
effective potential (see e.g. Ryan 1995) 



V(p)- 1 



£^^00 + 2ELg,^ + L^g„ 



(19) 



where, for circular orbits, the energy E and angular momen- 
tum L are determined by the conditions V = and dV/dp = 
(see Eqs. 22-23). The frequencies at the ISCO's location (de- 
termined by the additional condition d-y/dp^ = 0) are of par- 
ticular interest. Thus, before starting the discussion of the 
frequencies, it is important to explore the ISCO parametric 
dependence. We report here, as standard in the literature, the 
physical ISCO radius given by ^Jg^ evaluated at the root of 
Eq. (19) that gives the coordinate ISCO radius. In the upper 
panel of Fig. 1 we plotted contours of constant ISCO radii as 
a function of the dimensionless angular momentum parame- 
ter j - J/M^^ and the star quadrupole moment M2, for the 
PRS solution. The use of the dimensionless parameter j in 
the horizontal axis allows to, qualitatively, relate deviations 
of the contour lines from vertical lines to the influence of 
the quadrupole moment. We can see that the ISCO radius 
decreases for increasing j and decreasing M2. A quantitative 
measurement of this influence could be derived from the effec- 
tive slope of the contour lines. We are interested in the com- 
parison with the Kerr geometry, so in the lower panel, we plot- 
ted contours of constant ratio /"isco.PRs/nsco.Ken- as a function 
of j and the difference between the quadrupole moment of the 
PRS solution M2,prs and the Kerr quadrupole M2,Ken - -ma^, 
i.e. M2.PRS - M2,KeiT - A^2,PRS + ma' - mk, see Eq. (11). Devi- 
ations from the Kerr geometry are evident. Negative values of 
the angular momentum correspond to the radii of the counter- 
rotating orbits obtained here through the change g,^ — > -gt^ 
(see discussion below). 

We stress that the accuracy of the PRS solution for describ- 
ing the ISCO radius of realistic NSs was akeady shown to be 
higher with respect to other analytic models (see Pachon et al. 
2006, for details). In Table 1 we compare the ISCO radius 
for two rapidly rotating NS, models 20 and 26, of Table VI of 
Pappas and Apostolatos (2012) for the EoS L. The lowest mul- 
tipole moments of the analytic models are fixed to the numeri- 
cal values obtained by Pappas and Apostolatos (2012). In the 
case of the Kerr solution, only Mq and J can be fixed, while 
M2, and S3 have values that depend on Mq and J and there- 
fore cannot be fixed. For the PRS solution with i = 0, Mq, / 
and M2 can be fixed while S 3 remains induced by the lower 
moments. We present also the ISCO radius obtained by fixing 
Mo, J, M2, as well as 53 in the PRS analytic exact model. 

In Figs. 1-6, we have fixed as an example Mq - m - 
1.88Mq - 2.78 km, and i = 0. We recafl that the quadrupole 
moment in the geometric units used here (km"') is related 
to the one in COS units by M^"^^ = (10'^c2/G)Mf ° = 
1.35 X 10'*^(Mf7km^) g cm-, and the mass of the Sun is 
- 1 .477 km. The dimensionless angular momentum j is 
obtained from the CGS values of J and Mq as j = cJ/{GM~^). 

It is appropriate to compare the range of values of j = J/M^ 
and M2 used in Figs. 1-6 with typical values of a NS. For 
the used mass Mq - I.88M0, Morsinkand SteUa (1999) ob- 
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Fig. 1 . — Upper panel: Contours of constant ISCO radius as a function of the 
dimensionless angular momentum parameter j = J/M^ and the quadrupole 
moment M2 for the PRS solution, for a compact object with mass Mq = m = 
I.88M0 = 2.78 km. Contours are labeled by the corresponding the value 
of the ISCO radius in km. Negative values of j depict the counter-rotating 
case and negative values of the quadrupole moment M2 correspond to oblate 
configurations. The values of M2 are in the range < M2 < 20 km^^ that 
con'esponds in CGS units to < M2 < -2.7 X 10'*'* g cm^, which covers the 
typical range of fast rotating NSs. Lower panel: Contours of constant ratio 
nsco.PRs/nsco.Kcn as a function of j and the difference Mi.prs - Mi.Kerr- 
The quadrupole moment difference is comprised in the range -2.7 X 10^ < 
M2 < 6.8 X lO''^ gem-. 

tained a quadrupole moment M2 - -5.3 x 10"^^ g cm- = 3.93 
km-^, with the latter value in geometric units, for a NS of 
angular rotation frequency - 290 Hz (rotation period 
of 3.45 milliseconds), corresponding to a dimensionless an- 
gular momentum j - J/Mq - 0.19, for the EoS L. For 
a fixed mass the quadrupole moment is an increasing func- 



M20 
M26 



19.87 



13.39 
17.16 



TSIT 
15.94 



19.28 
19.65 



18.99 
19.54 



TABLE 1 

Comparison of the ISCO radius for the selected NS models 20 and 26 of 
Table VI of Pappas and Apostolatos (2012) for the EoS L. Model 20: 
Mo = 4.167 KM (2.82Mo), j = J/M^ = 0.70, M2 = -79.8 km' (-1.08 x lO'*' 

G CM^) AND S3 = -401.0 km". Model 26: Mq = 4.36 km (2.95Mo), 
j = JImI = 0.56, M2 = -45.2 km' (-6.10 x lO""* g cm^) and 53 = -170.0 

km"*. The subscript N stands for the numerical calculation of 
Pappas and Apostolatos (2012) and SS stands for the Shibata and Sasaki 
(1998) approximated ISCO radius expression. 



tion of i because an increasing of the angular momentum at 
fixed mass results in an increasing of the oblateness (eccen- 
tricity) of the star, and so the quadrupole moment. Based 
on this fact, it is clear that not all the values of the {M2,j) 
pairs of quadrupole and angular momentum depicted in, e.g.. 
Fig. 1 are physically meaningful. The maximum rotation 
rate of a neutron star taking into account both the effects 
of general relativity and deformations has been found to be 
Vs,max = lO45(Mo/M0)'/2(lOkm/7;)^/- Hz, largely indepen- 
dent on the EoS (see Lattimer and Prakash 2004, for details). 
Corresponding to this maximum rotation rate, the angular mo- 
mentum is /max - 2;7rVs,niax^ ~ 6.56 X 10'*^/45 g cm" s"', 

and y„ax = GJ^,^/(cMI) ~ O.lAhiUMolM^f , where hi 
is the moment of inertia of the NS in units of 10"^^ g cm^. 
The fastest observed pulsar is PSR J1748-2246ad with a ro- 
tation frequency of 716 Hz (Hessels et al. 2006), which con- 
strains the mass of the NS to Mq > O.47(/?/lOkm)^M0, and 
j ~ 0.51I45/(Mo/Mq)^, which becomes ; ~ O.26/45 for a 
canonical NS of Mq = 1 .4Mq. 

3.1. Keplerian Frequency 

Now we turn into the frequencies analysis. For stationary 
axially symmetric spacetimes, the frequency of Keplerian or- 
bits is given by (see e.g. Ryan 1995) 



(20) 



where a colon stands for partial derivative with respect to the 
indicated coordinate and and '-' stands for corotating and 
counter-rotating orbits, respectively. 

For the case of static spacetimes, i.e. for co = and there- 
fore g,^ = 0, Qk = ± ^J-84,4,,p8ii,p/g4,4,.p and the energy E and 
angular momentum L per mass /i of the test particle can be 
expressed in terms of the metric tensor components (see e.g. 
Ryan 1995), 



E 



L 



(21) 



From here, it is clear that taking the negative branch of the 
root for in Eq. (20) is equivalent to studying a particle 
with opposite angular momentum, i.e. Lcoum-rot - -i'co-rot- 
Thus, in the static case the magnitude of the energy and angu- 
lar momentum are invariant under the change ^ -^k- 

Now we consider the case of stationary space times, w ?i 0. 
The energy E and angular momentum L per mass fi are, in this 
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case, given by (see e.g. Ryan 1995) 

E -gl, - gl4>^K 



L _ g,^ + g00^K 



(22) 
(23) 



The counter-rotating condition expressed by the negative 
branch of Eq. (20), can be generated by the change g,^ 
-gtip, which seems to be a more physical and transparent con- 
dition. In contrast to the static case, the counter-rotating or- 
bit has now different energy and different magnitude of the 
angular momentum due the presence of the dragging of in- 
ertial frames, characterized by the metric component g,^ (cf. 
Eq. (26) below). In a nutshell, the dynamics of counter- 
rotating orbits of a test-particle can be derived, starting from 
the positive branch of Eq. (20), by considering a spacetime 
with g,^ -g,^. 

For the vacuum case, a similar analysis as the one develo- 
ped by Herrera et al. (2006), clearly shows that the change in 
the global sign of g,^ is achieved by changing not only the 
angular momentum of the star, J -J, but all the rota- 
tional multipolar moments. For the Kerr metric this change 
is obtained by changing the sign of the parameter a (see Ap- 
pendix B) while in the PRS solution we need additionally 
change the sign of the parameter s associated to differential 
rotation, i.e., by changing a ^ -a and s -s^. 

Once we have clarified this important issue about the co- 
rotating and counter-rotating orbits, we proceed to analyze 
the functional dependence of the Keplerian frequency on the 
multipole moments. In the upper panel of Fig. 2 we plotted 
contours of constant Keplerian frequency for the PRS solu- 
tion, vk.prs - i^K,PRs/(27r), as a function of the dimensionless 
angular momentum parameter j and the quadrupole moment 
Mo.pRS, at the ISCO radius. It can be seen that the influence of 
the quadrupole moment is non-negligible, as evidenced from 
the departure of the contour lines from vertical lines. The Ke- 
plerian frequency grows for increasing / and M2. In the lower 
panel, we plotted contours of constant ratio vk.prs / vk.kbh as a 
function of j and the difference between the quadrupole mo- 
ment of the PRS solution, M2,prs, and the Kerr quadrupole, 

M2,KeiT- 

It is appropriate to recall here that because the Keplerian 
as well as the other frequencies calculated below are evalu- 
ated using formulas in the coordinate frame, see for instance 
Eq. (20), they must be evaluated at coordinate radii p and 
not at physical radii given by yg^^. In the specific case of 
the ISCO the frequencies are evaluated at the radius that si- 
multaneously solves the equations V - Q, dV/dp = 0, and 
d^V/dp^ - 0, where V is the effective potential (19). 

3.2. Oscillation and Precession Frequencies 

The radial and vertical oscillation (or epicyclic) frequencies 
are the frequencies at which the periastron and orbital plane 
of a circular orbit oscillates if we apply slightly radial and 
vertical perturbations to it, respectively. According to Ryan 
(1995), in stationary axially symmetric vacuum spacetimes 
described by the Weyl-Papapetrou metric (1), the radial and 

^ For the vacuum case, in the solution by Manko et al. (2000), the sign 
change of g,^ is obtained after performing simultaneously the replacements 
a —> -a and b —> -b. 
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Fig. 2. — Upper panel: Contours of constant vk (in kHz) as a function of the 
the dimensionless angular momentum parameter j = J/M^ and quadruple 
moment M2 for the PRS solution, at the ISCO radius, for a compact object 
with mass Mq = in = I.88M0 = 2.78 km. Lower panel: Contours of constant 
ratio Vk prs / Keir as a function of ; and the difference M2 prs - M2 Kcm at 
the ISCO radius. 



vertical epicyclic frequencies can be obtained as 



..OS 



1 

2^ 



{g„ + gt^siy I -2- 



2(g„ -1- g,^n){g,^ + g^^Q.) I — 

\P /,aa 



+ {g,^+g^^nf[^^ 



1/2 



(24) 
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and the corresponding periastron (v^) and nodal (vf) preces- 
sion frequencies as 



Vk 



..OS 



(25) 



where a - {p,z}, respectively, and vk - Q.k/{27t) is the Keple- 
rian orbital frequency with Qk given by Eq. (20). 

In the upper panel of Fig. 3, we plotted contours of con- 
stant nodal precession frequency vf at the ISCO radius as 
a function of j - J/M^ and M2 for the PRS solution, at 
the ISCO radius. We can see now that the influence of the 
quadrupole moment is quite important. The nodal precession 
frequency increases for increasing J and decreasing M2, at 
fixed Mq. In the lower panel we plotted contours of constant 
ratio T^"^ /^"^ Ken-' '■^^ ISCO radius, as a function of j and 
the difference Mz.prs - A^2,KeiT, in order to evidentiate devia- 
tions from the Kerr solution. The radial oscillation frequency 
vanishes at the ISCO radius and therefore at such loca- 
tion the radial precession frequency equals the Keplerian fre- 
quency, whose contours have been plotted in Fig. 2. 

In Figs. 4 and 5 we plotted the nodal precession frequency 
and the radial oscillation frequency as a function of 
the Keplerian frequency vk, respectively, for both the Kerr 
and PRS solutions. As an example, we have shown the 
results for rotating NS models 20 and 26 of Table VI of 
Pappas and Apostolatos (2012), for the EoS L. The lowest 
multipole moments of the PRS solution Mq, J, M2, and 
53 have been fixed to the numerical values obtained by 
Pappas and Apostolatos (2012). In the case of the Kerr so- 
lution, only Mo and J can be fixed, while M2, and 53 have 
values induced by the lower moments Mq and J. For the PRS 
solution with i = 0, Mq, J and M2 can be fixed while S 3 can- 
not be fixed and depends on the lower moments. The results 
for the PRS analytic model obtained by fixing Mq, J, M2, as 
well as S 3 are also shown. 

The deviations of the quadrupole and current octupole mo- 
ments given by the Kerr solution from the numerical values 
of Pappas and Apostolatos (2012) can be used to show the 
low accuracy of the Kerr solution to describe fast rotating NSs. 
The accuracy of the PRS solution in describing the ISCO radii 
of these two models has been shown in Table 1 of Section 3. 

In Figs. 4 and 5 we can see the differences of the v^-vk 
and Vp ^-vk relations between the Kerr and PRS solutions for 
realistic NS models. The deviations of the Kerr solution, es- 
pecially at fast rotation rates, are evident because of the influ- 
ence of the deformation (quadrupole M2) of the star as well 
as, although in less proportion, of the octupole current 5 3. In 
general, we observe that the larger the angular momentum, the 
poorer the performance of the predictions of Kerr solution. 

We have also shown in Figs. 4-5 the influence of the current 
octupole 5 3 in the determination of the precession and oscil- 
lation frequencies. We found that the eff'ect of 53 is only ap- 
preciable for the fastest models. The minor influence, in this 
case, of the current octupole 5 3 is expected from the small 
values of the parameter s needed to fit the numerical values of 
Pappas and Apostolatos (2012). Clearly, larger values of the 
parameter s needed to fit realistic values of 5 3 will enhance 
as well deviations from the Kerr spacetime. 

The effects of a multipolar structure that deviates from 
the one of the Kerr geometry on the various quantities 
analyzed here are relevant for instance in the RPM of 
the QPOs observed in LMXBs (see e.g. Stella and Vietri 
(1998); Morsink and Stella (1999); Stella ef a/. (1999); 
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Fig. 3. — Upper panel: v? (in kHz) as a function of the the dimensionless 
angular momentum parameter ; = J/M^ and quadruple moment M2 for the 
PRS solution, at the ISCO radius, for a compact object with mass Mq = m = 
I.88M0 = 2.78 km. Lower panel: Contours of constant ratio i'^prs /"-Ken 
a function of j and the difference Mi.prs - M2_Kcrr. at the ISCO radius. 

Stella and Vietri (1999) and Section 6, for details). 

3.3. Dragging oflnertial Frames 

It is known that a prediction of general relativity is that a ro- 
tating object makes a zero angular momentum test particle to 
orbit around it, namely it drags the particle into the direction 
of its rotation angular velocity; such an effect is called drag- 
ging of inertial frames or Lense-Thirring effect. Consequently, 
oblique particle orbit planes with respect to the source equa- 
torial plane will precess around the rotation axis of the object. 
In stationary axially symmetric spacetimes described by the 
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Fig. 4. — Nodal precession frequency y? versus Keplerian frequency vk 
given by the Kerr and PRS analytic solutions. The lowest multipole mo- 
ments have been fixed from the rotating NS models 20 (red curves) and 26 
(blue curves) of the Table VI of Pappas and Apostolatos (2012) for the EoS 
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Fig. 5. — Radial oscillation frequency v" versus Keplerian frequency vk 
given by the Kerr and PRS analytic solutions. The lowest multipole mo- 
ments has been fixed from the rotating NS models 20 (red curves) and 26 
(blue curves) of the Table VI of Pappas and Apostolatos (2012) for the EoS 
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metric (1) the frame dragging precession frequency is given 
by (see e.g. Ryan 1995) 



Vlt = - 



7.n I 



(26) 



Many efforts have been done to test the predictions of 
general relativity around the Earth such as the analysis of 
the periastron precession of the orbits of the LAser GEO- 
dynamics Satellites, LAGEOS and LAGEOS II, (see e.g. 
Lucchesi and Peron 2010) and the relativistic precession of 
the gyroscopes on-board the Gravity Probe B satellite (see 
Everitt ef flZ. 2011, for details). The latter experiment mea- 
sured a frame dragging effect within an accuracy of 19% with 
respect to the prediction of general relativity. 

The smallness of this effect around the Earth makes such 



measurements quite difficult and has represented a multi year 
challenge for Astronomy. The frame dragging precession in- 
creases with the increasing of the angular momentum of the 
rotating object and therefore a major hypothetical arena for 
the searching of more appreciable Lense-Thirring precession 
is the spacetime around compact objects such as BHs and NSs. 
The much stronger gravitational field of these objects with re- 
spect to the Earth one allow them to attain much faster angular 
rotation rates and so larger angular momentum. 

Stella and Vietri (1998) showed how, in the weak field slow 
rotation regime, the vertical precession frequency vF (orbital 
plane precession frequency) can be divided into one contribu- 
tion due to the Lense-Thirring precession and another one due 
to the deformation (non-zero quadrupole moment) of the ro- 
tating object, both of them comparable from the quantitative 
point of view. These frequencies could be in principle related 
to the motion of the matter in the accretion disks around BHs 
and NSs and thus particularly applicable to LMXBs. For fast 
rotating NSs and BHs the frequency at which the orbital plane, 
and so the frame dragging precession frequency, can reach val- 
ues of the order of tens of Hz (see e.g. Stella and Vietri (1998) 
and Figs. 3 and 4). 

Thus, it is clear that an observational confirmation of the 
relativistic precession of matter around either a NS or a BH 
will lead to an outstanding test of the general relativity in the 
strong field regime and, at the same time, an indirect check 
of the large effects of the frame dragging in the exterior space- 
time of compact objects (see e.g. Morsink and Stella 1999, for 
details). 

Although making independent measurements of the frame 
dragging effect around BHs and NSs is a very complicate task, 
it is important to know the numerical values of the precession 
frequency due to the frame dragging with respect to other rel- 
ativistic precession effects, e.g. geodetic precession. In addi- 
tion, it is important to know the sensitivity of the precession 
frequency to the object parameters such as mass, angular mo- 
mentum, quadrupole, and octupole moment. 

In the upper panel of Fig. 6 we plotted contours of con- 
stant frame dragging frequency vlt for the PRS solution, at 
the ISCO radius, as a function of the the angular momentum 
per unit mass //Mo and the quadruple moment M2, for a com- 
pact object mass Mo = m - I.88M0. Correspondingly, in 
the lower panel of Fig. 6, we show the differences between 
the frame dragging precession frequency as predicted by the 
Kerr and PRS solutions, at the ISCO radius, as a function 
of j - JIM^ and the difference between the quadrupole mo- 
ments, Mi.pRS - M2,Ken- 

The influence of the quadrupole moment in the determina- 
tion of the frame dragging frequency is evident; the frequency 
Vlt given by a NS is generally smaller than the one given by a 
BH as can be seen from the value of the ratio vlt.prs /TLX.Ken- < 
1 obtained for configurations with a quadrupole moment that 
deviates with respect to the one given by the Kerr solution, 
namely for Mt prs - Mt kmt = M2 prs + ma^ - mk + 0, see 
Eq.(ll). 

It is also worth mentioning that frame dragging precession 
can be affected as well by the presence of electromagnetic 
fields (Herrera et al. 2006, see) and further research in this 
respect deserves the due attention. 

4. ACCURACY OF RYAN'S ANALYTIC FORMULAS 

Following a series expansion procedure in powers of 1/p, 
Ryan (1995) found that the periastron (radial) and nodal (ver- 
tical) precession frequencies, Vp and vf given by Eq. (24), can 
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Fig. 6. — Upper panel: Contours of constant vlt (in Hz) as a function of the 
the angular momentum per unit mass j = J/M^ and the quadruple moment 
M2 for the PRS solution, at the ISCO radius, for a compact object with mass 
Mo = m = I.88M0 = 2.78 km. Lower panel: Contours of constant ratio 
Vlt PRS / I'LT Kerr as a function of jO and the difference M2 prs - M2 Kerr . at the 
ISCO radius. 



be written as a function of the Keplerian frequency vk as 
VK Ml [2 2 Ml) Ml 



(21 
J 



5? 21 Mj 



' M* 



^6 + 



405 ^ 2243 S] 



84 Ml 



-48^-: 
Ml 

661 M2 
' 14 Ml 



. S1M2 
'"0 



21 



^ 15 M4' 



(27) 



and 



Vk Ml 2 Ml 



7-^+3 



M2 

~3 



11 



Ml "m: 



39Ml _ 15 M4 

Q A ^Jr5 



^7 ^ 



\_ 153 M2 

28" m4 



153 52 



8 Ml 



4 



(28) 



where 'V = (2;rMoVK)'^^ \MqMi,Ma\ are the lowest three 
mass moments and, [5 1,53], are the lowest two current mo- 
ments. For the PRS solution in the vacuum case, M4 - 
m{a' — licr + ]c + 2as). 

The above formulas are approximate expressions of the pe- 
riastron and nodal precession frequencies in the weak field 
(large distances from the source) and slow rotation regimes. 
We should therefore expect that they become less accurate at 
distances close to the central object, e.g. at the ISCO radius, 
and for fast rotating objects. However, such formulas are an 
important tool to understand the role of the lowest multipole 
moments on the values of the relativistic precession frequen- 
cies, such as the importance of the higher multipole moments 
at short distances and high frequencies as can be seen from 
Eqs. (27-28). 

At high frequencies, for instance of the order of kHz, devi- 
ations from the above scaling laws are appreciable. In Figs. 7 
and 8 we compare the radial precession and vertical oscilla- 
tion frequencies, v^ and v?^, as a function of the Keplerian fre- 
quency Vk, as given by the full expressions (24) for the PRS 
solution and by the approximate formulas (27) and (28), re- 
spectively.^ The lowest multipole moments Mq, 7, M2, and 5 3 
of the PRS solution have been fixed to the values of two mod- 
els of the Table VI of Pappas and Apostolatos (2012); Model 
2 with Mo = 2.071 km (I.4O2M0), ; = 0.194, M2 = -2.76 
km^^ (3.73 X 10'*^ g cm^), 53 = -2.28 km^ and Model 20 with 
Mo = 4.167 km (2.82M0), ; = J I Ml = 0.70, M2 = -79.8 
kmV-l-08x gcm^)wASi = -401.0 km^. 

In the v?^-vk relation, the blue dotted curve depicts the 
contribution from the angular momentum (we plot the series 
(28) up to 'V^), for the blue dot-dashed curve we added the 
first contribution from the quadrupole moment M2 (we cut 
the series at 'V'*), for the dashed blue line we added the first 
contribution from the octupole mass-current (series expansion 
up to 'V'') and finally in the continuos blue line we consider 



and V? frequencies are very similar, we de- 



Because the scale of the rp 
cided to plot in Fig. 7 and v?^ 
clear comparison with the PRS solution in a single figure. 



whose scales are different allowing a more 
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Fig. 7. — Comparison of the vP^-vk and v^-vk relations given by the PRS 
solution and the approximate expressions (27-28) derived by Ryan (1995). 
The lowest multipole moments Mo, J, M2, and S3 have been tixed to the 
values of the Model 2 of the Table VI of Pappas and Apostolatos (2012): 
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Fig. 8. — Comparison of the v9^-vk and v^-vk relations as given by the PRS 
solution and the approximate expressions (27-28) derived by Ryan (1995). 
The lowest multipole moments Mo, ^2, and S3 have been fixed to the 
values of the Model 20 of the Table VI of Pappas and Apostolatos (2012): 
Mo = 4.167 km (2.82Mo), i = 7/M^ = 0.70, M2 = -79.8 km^ (-1.08x10''' 
g cm^) and S3 = -401.0 km'' 

contributions for higher multipole moments and stop the se- 
ries at the order 'V^, not shown in Eqs. (28). For this case, we 
can see that Ryan's expressions clearly tend, from the bottom, 
to the exact result (continuous black curve) obtained by using 
the PRS solution. 

For the analysis of the Vp-VK relation we followed the same 
procedure as described above. In this case, the Ryan's ex- 
pressions tend from the top to the exact result, the continuous 
black curve, represented by the PRS solution. It is interesting 
to see that the introduction of the octupole moment (dashed 
red line) makes the approximation to deviate from the exact 
result, however by including more terms the accuracy is en- 
hanced. As can be seen from Figs. 7 and 8 the quantitative 
accuracy of the Ryan's approximate formulas in the perias- 
tron precession frequency is less than the one obtained in 
the vertical oscillation frequency . 

The importance of the high-order multipole moments such 
as the quadrupole and the octupole moments is evident in the 
high-frequency regime. This is in line with the results shown 



in Figs. 2-3 and in Figs. 4—5. We can see from Figs. 7 and 
8 that the Ryan's approximate formulas describe more accu- 
rately the Model 2 than the Model 20. The reason is that, as 
we mentioned above, we should expect a better accuracy of 
the series expansions from low to moderate moderate rotation 
rates and consequently the same occur for the quadrupole de- 
formations. It is clear that there are appreciable differences 
both in rotation and deformation between the two selected 
models; we recall also that the rotation frequency of the star 
can be expressed as a function of the dimensionless j parame- 
ter as V, ^GjMllilncI) = 1.4(M/MO)-//45 kHz. 

It is noteworthy that we have checked that the Ryan's series 
expansions, Eqs. (27) and (28), fit quite accurately the exact 
results if taken up to order 'V'". In particular the values of 
the vertical oscillation and precession frequencies are fit better 
than the corresponding radial ones. For the Model 2 the radial 
oscillation frequency is well fitted by the Ryan's expression 
up to Keplerian frequencies of order ~ 1 .2 kHz while, for the 
Model 20, the approximate formulas break down at a lower 
value ~ 0.7 kHz. These results are of particular relevance be- 
cause it makes possible the extraction of the object parameters 
(for instance the lowest multipoles up to S 3) by the fitting of 
the observed QPO frequencies in LMXBs, providing they are 
indeed related to the precession and oscillation frequencies of 
matter in the accretion disk (see Section 6, for details) and for 
Keplerian motion not exceeding a few kHz of frequency. 

5. ACCURACY OF PRS SOLUTION 

We turn now to analyze the behavior of the Kerr and 
PRS solutions in predicting results for the Keplerian, frame 
dragging, and vertical oscillation frequencies, for realistic 
NSs. In particular, we compare their predictions with the 
frequencies calculated by Morsink and Stella (1999). Since 
Morsinkand Stella (1999) did not include the values of the 
octupole current moment 5 3, here we set i = in Eq. (10) 
for the PRS solution. For the sake of comparison, we choose 
the results derived by Morsink and Stella (1999) for the EoS 
L, because for this EoS the highest rotating parameter j and 
quadrupole moment M2 were found. In addition, the stiff- 
ness of such an EoS allows the maximum mass of the NS 
to be larger than the highest observed NS mass, Mq - 1 .97 ± 
O.O4M0, corresponding to the 317 Hz (3.15 milliseconds rota- 
tion period) pulsar J16I4-2230 (see Demorest et al. 2010, for 
details). 

This regime of high j and M2 in realistic models is partic- 
ularly interesting to test the deviations of the Kerr solution in 
the description of NS signatures as well as to explore the ac- 
curacy of the PRS solution. In Table 2, we present the results 
for four different sets of the star spin frequency v,, namely 
V, = 290 Hz (Ml and M2), v, = 360 Hz (M3 and M4), 
V, = 580 Hz (M5 and M6) and v, = 720 Hz (M7 and M8). 

In Table 2, we clearly observe that the results predicted 
by the PRSj=o solution for the Keplerian and frame-dragging 
frequencies are in excellent agreement with those calculated 
by Morsink and Stella (1999) for even highly massive, ro- 
tating and deformed models such as the model M7 with 
Mo = 2.17Mo, = 0.51 and M2 = -39.420- We notice that 
Morsink and Stella (1999) reported some configurations with 
negative values of (see Table 2). We advance the possibility 
that this is due to instabilities of the numerical code that occur 
when the ISCO radius is located very close or inside the sur- 
face of the object. Thus, the values of the frequencies given by 
the analytic solution in these cases are to be considered predic- 
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0.19 
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15.4 


14.90 


15.42 


1.31 


1.363 


1.304 


39.7 


42.248 


38.476 


18.6 


39.697 


22.040 


M2 


2.71 


0.14 


-3.0 


22.2 


22.16 


22.23 


0.90 


0.906 


0.902 


19.6 


19.676 


19.493 


17.2 


18.809 


17.629 


M3 


1.94 


0.24 


-8.2 


15.6 


14.89 


15.63 


1.29 


1.380 


1.296 


49.8 


57.001 


49.804 


17.5 


52.608 


26.197 


M4 


2.71 


0.18 


-4.8 


21.8 


21.62 


21.74 


0.93 


0.937 


0.931 


26.1 


27.245 


26.833 


21.9 


25.670 


23.635 


MS 


2.07 


0.40 


-23.1 


16.3 


14.18 


16.06 


1.26 


1.514 


1.289 


84.3 


125.75 


88.905 


-10.5 


109.04 


31.140 


M6 


2.72 


0.30 


-13.9 


20.6 


20.05 


20.45 


1.01 


1.041 


1.015 


53.5 


57.467 


54.391 


35.8 


51.861 


42.854 


M7 


2.17 


0.51 


-39.4 


17.0 


13.58 


16.53 


1.22 


1.637 


1.269 


106.7 


201.52 


116.25 


-51.8 


166.55 


29.804 


M8 


2.73 


0.39 


-24.4 


19.8 


18.85 


19.65 


1.06 


1.136 


1.077 


78.8 


90.821 


80.895 


38.4 


79.079 


57.187 



TABLE 2 

ISCO RADIUS r+, KePLERIAN frequency Vk, frame-dragging (LeNSE-ThIRRING) frequency VlT, and vertical precession frequency V? OF THE CO-ROTATING ORBITS 
CALCULATED NUMERICALLY BY MoRSINK AND StELLA (1999) (uPPER INDEX MS) AND COMPARISON WITH THE CORRESPONDING PREDICTED VALUES GIVEN BY THE KeRR (uPPER 

INDEX Kerr) and the PRS.5=o solution (upper index PRS). The quadrupole moment M2 have been normalized for convenience to the value Qo = 10"*' g cm^. 



tions to be tested for future numerical computations. This fact 
can be checked within the calculations of Morsink and Stella 
(1999) by exploring the properties of counter-rotating orbits 
which produce in general ISCO radii larger than the ones of 
the corotating ones. In Table 3, we depicted the results in the 
counter-rotating case where we can notice an improvement 
of the accuracy of the PRS solution with respect to the co- 
rotating case. 

In this line, we consider worth performing numerical com- 
putations of the precession and oscillation frequencies of par- 
ticles around realistic NSs in a wider space of parameters and 
using up-to-date numerical techniques which will certainly 
help to establish and elucidate more clearly the accuracy of 
analytic models. It is also appropriate recalling the recent re- 
sults of Pappas and Apostolatos (2012) on the computation of 
the general relativistic multipole moments in axially symmet- 
ric spacetimes. 

6. THE RELATIVISTIC PRECISION MODEL 

The X-ray light curves of LMXBs show a variability from 
which a wide variety of QPOs have been measured, expand- 
ing from relatively low ~ Hz frequencies all the way up to 
high ~ kHz frequencies (see e.g. vanderKlis 1995, for de- 
tails). In particular, such frequencies usually come in pairs 
(often called twin peaks), the lower and upper frequencies, v; 
and Vh respectively. BHs and NSs with similar masses can 
show similar signatures and therefore the identification of the 
compact object in a LMXB is not a simple task. If the QPO 
phenomena observed in these systems are indeed due to rela- 
tivistic motion of accretion disk matter, the knowledge of the 
specific behavior of the particle frequencies (e.g. rotation, os- 
cillation, precession) in the exterior geometry of NSs and BHs 
becomes essential as a tool for the identification of the nature 
of the compact object harbored by a LMXB. 

It is not the scope of this work to test a particular model 
for the QPO phenomenon in LMXBs but instead to show the 
influence of the high multipole moments on the orbital mo- 
tion of test particles especially the role of the quadrupole mo- 
ment which is of particular interest to differentiate a NS from a 
BH. There are in the literature several models that describe the 
QPOs in LMXBs through the frequencies of particles around 
the compact object, and for a recent review and comparison of 
the different models we refer to the recent work of Lin et al. 
(2011). In order to show here the main features and differ- 
ences between the Kerr and the PRS solutions we shall use 
the Relativistic Precession Model (RPM). 

The RPM model identifies the lower and higher (often 
called twin-peaks) kHz QPO frequencies, v/ and v;,, with 
the periastron precession and Keplerian frequencies, namely 
V/ = Vp and Vh = vk, respectively. The so-called horizon- 
tal branch oscillations (HBOs), which belong to the low fre- 
quency QPOs observed in high luminosity Z-sources (see e.g. 



vanderKlis 1995, for details), are related within the RPM 
model to the nodal precession frequency v? of the same orbits 
(Morsink and Stella 1999, see). We will use here in particular 
the realistic NS models of Morsink and Stefla (1999) for the 
EoS L. 

One of the salient features of the RPM model is that in the 
case of the HBO frequencies, the relations inferred from the 
first term of the expansions (27) and (28) 

yK = 3-^/5 (2;r)-2/5,„-2/5(^)3/5^ (29) 
yP = (2/3)«/V/5jm'/5(,^)6/5, (30) 

which implies a nodal precession frequency proportional to 
the square of the Keplerian frequency has been observed in 
some sources, for instance in the LMXB 4U 1728-34 (see 
Ford and van der Klis 1998, for details). In addition, 6/5 
power law relating the nodal and periastron precession fre- 
quencies can explain (see Stella ef a/. 1999) the correlation 
between two of the observed QPO frequencies found in the 
fluxes of NSs and BHs LMXBs (see Psaltis et al. 1999, for 
details). This fact provides, at the same time, a significant test 
of the Ryan's analytic expressions. 

It is interesting to analyze the level of predictability of the 
precession and oscillation frequencies on particular astrophys- 
ical sources. In Fig. 9 we show the v/-v;, relation within the 
RPM model, namely versus vk for the models M1-M8 of 
Table 2. In the upper panel we show the results for the PRS 
solution while, in the lower panel, we present the results for 
the Kerr solution. We have indicated the QPO frequencies ob- 
served in the sources GX 5-1 (see e.g. Wijnands et al. 1998; 
Jonkeref aZ. 2002), 4U 1735-44 (see e.g. Foid etal. 1998), 
4U 1636-53 (see e.g. Wijnands ef a/. 1997), Sco XI (see 
e.g. van der Klis et al. 1996), GX 17-2 (see e.g. Homan et al. 
2002), GX 340+0 (see e.g. Jonker et al. 2000), Cir XI (see e.g. 
van der Klis ef fl/. 1996), 4U 0614-H091 (see e.g. Ford etal. 
1997), and 4U 1728-34 (see e.g. Strohmayer ef al. 1996). 

Both the upper and lower panels of Fig. 9 have been plot- 
ted using the same frequency scales in order to aid the iden- 
tification of the differences between the Kerr and the PRS 
solutions. One can notice that all the solid curves in the 
Kerr solution (lower panel of Fig. 9) are outside the range of 
the observed QPO frequencies exemplified, while all dashed 
and solid curves of the PRS are inside the QPO range. It 
is then clear that making a fit of the observed QPO frequen- 
cies of the selected LMXBs of Fig. 9 will necessarily require 
a different choice of parameters in the Kerr and PRS solu- 
tions. Therefore, conclusions for instance on the NS parame- 
ters (e.g. mass, angular momentum, quadrupole deformation) 
based on fitting QPOs using the Kerr geometry will deviate 
from the actual parameters (see e.g. Laarakkers and Poisson 
1999, for details), extractable more reliably from a more com- 
plex geometry, such as the PRS one, that allows a better es- 
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Model 


r"S[km] 




^PRS 


v^^[kHz] 


K 




'1 


^Kcrr 
I T 


yPRS 
IT 






^PRS 


Ml 


18.8 


18.35 


18.73 


-V59— 


1.023 


0.997 


21.7 


22.61 


21.39 


29.0 


23.88 


29.70 


M2 


25.9 


25.82 


25.88 


0.73 


0.734 


0.732 


12.4 


12.43 


12.36 


13.6 


12.95 


13.46 


M3 


19.9 


19.39 


19.89 


0.93 


0.960 


0.928 


24.0 


25.79 


24.04 


32.9 


27.61 


34.20 


M4 


26.4 


26.33 


26.42 


0.71 


0.715 


0.712 


14.7 


15.08 


14.94 


16.4 


15.88 


16.57 


M5 


23.3 


22.22 


23.24 


0.77 


0.816 


0.768 


28.2 


32.59 


28.78 


41.3 


36.31 


44.65 


M6 


28.2 


27.94 


28.18 


0.65 


0.660 


0.652 


20.5 


21.19 


20.70 


24.1 


23.03 


24.51 


M7 


25.9 


24.37 


25.78 


0.67 


0.731 


0.678 


28.9 


34.62 


29.64 


43.4 


39.56 


47.90 


M8 


29.7 


29.18 


29.58 


0.61 


0.623 


0.611 


23.2 


24.40 


23.48 


28.4 


27.12 


29.22 



TABLE 3 

Same as in Table 2, but for the counter-rotating case. 



0.5 



A 


GX 5-1 


+ Sco XI 


■ Cir XI 


T 


4U 1735 -44 


• GX 17-2 


X 4U 0614+091 


♦ 


4U 1636-53 


GX 340 +0 


• 4U 1728-34 




Fig. 9. — Periastron oscillation frequency, y^^, as a function of the Keplerian 
frequency vk for the NS realistic models in Table 2. We indicate the QPO 
frequencies observed in the sources GX 5-1, 4U 1735^4, 4U 1636-53, Sco 
XI, GX 17-2, GX 340+0, Cir XI, 4U 0614+091, and 4U 1728-34. The 
solid curves depict the results for the models Ml (solid) and M2 (dashed) 
with red lines, for the models M3 (solid) and M4 (dashed) with blue lines, for 
the models M5 (solid) and M6 (dashed) with green Unes while orange Unes 
stands for the results from models M7 (soUd) and M8 (dashed). In the upper 
panel we present the results derived from the PRSj^o solution while in the 
lower panel we present the results for the Kerr solution. In the lower panel 
we have added, to guide the eye, the inner red dashed and outer red sohd 
curves of the upper panel using black lines. 

timate for instance of the quadmpole moment of a compact 
star. 

In Fig. 9 we show the relation v? versus vk for the mod- 
els M1-M8 of Table 2. For the sake of comparison we show 
the low frequency branch observed in the LMXB 4U 1728-34 
(see Ford and van der Klis 1998, for details). From the analy- 
sis of the pulsating X-ray flux it turns out that very likely the 
spin frequency of the NS in 4U 1728-34 is ~ 363 Hz (see 
Strohmayer ef fl/. 1996, for details). Thus, the models M3 
(Mo = 1.94Mo,; = 0.24) and M4 (Mq = 2.71Ma,; = 0.18)in 
Table 2 that correspond to a NS of spin frequency 360 Hz are 
of particular interest for the analysis of this source. It was sug- 
gested by Stella et al. (1999); Stella and Vietri (1999) that the 
low frequency observed in 4U 1728-34 are Ukely to be due 




Fig. 10. — Nodal precession frequency, y. , as a function of the Keplerian 
frequency vk for the NS realistic models in Table 2. The convention is as 
Fig. 9. We indicate the QPO frequencies observed in the LMXB 4U 1728-34 
(see Ford and van der Klis 1998). The black curves indicate the 2v?-vk rela- 
tion for the models M3 and M4 (solid and dashed) following the suggestion 
of Stella et al. (1999); Stella and Vietri (1999). 

to excitations of the second harmonic of the vertical motion 
and therefore a better fit of the lower-higher QPO frequencies 
of 4U 1728-34 (and of similar sources) will be obtained for 
the relation 2v?-vk. The black curves in Fig. 10 indicate the 
2vf-VK relation for the models M3 and M4 (solid and dashed) 
following the above suggestion. Although the improvement 
of the fit is evident, we notice that the NS parameters that 
correctly reproduce the features of 4U 1728-34 are likely in 
between the models M3 and M4. 

7. CONCLUDING REMARKS 

We have done an extensive comparison of the orbital mo- 
tion of neutral test particles in the PRS and Kerr spacetime 
geometries. In particular we have emphasized on the Keple- 
rian and frame-dragging frequencies, as well as the precession 
and oscillation frequencies of the radial and vertical motions. 

We have evidentiated the differences in this respect between 
the Kerr and PRS solution, especially in the rapid ~kHz ro- 
tation regime. Such differences are the manifestation of the 
influence of the high order multipole moments such as the 
quadrupole and octupole. 
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The analysis of the deviations between the Kerr and PRS 
features for given mass and angular momentum of a source 
studied in this work are useful to distinguish the signatures 
between BHs and NSs, which relevant to establish a separa- 
trix for the identification of the compact objects harbored in 
X-Ray Binaries. In the case of BH candidates, these results 
might become important for testing the no-hair theorem of 
BHs (see e.g. Johannsen and Psaltis 201 1). Equally important, 
the application of the precession and oscillation frequencies 
to the explanation of QPOs in LMXBs possessing a NS, can 
unveil information on the NS parameters, leading to a possi- 
ble identification of the behavior of the nuclear matter EoS at 
supranuclear densities. In this line, the identification of the ro- 
tation frequency of NSs in LMXBs from the pulsating X-ray 
flux Vburst, e.g. the case of 4U 1728-34 (Ford and van der Klis 
1998), 4U 1916-053 (Galloway et al. 2001) and more recently 
the case of IGR J17191-2821 (Altamirano et al. 2010), will 
certainly help to constrain QPO models as well as the NS 
parameters. Additional information coming from the photo- 
spheric radius expansion phenomena observed in these sys- 
tems (see e.g. Muno et al. 2001, for details) during their tran- 
sient activity with Super-Eddington emission can become of 
paramount importance if combined with the QPO informa- 
tion. 

The generalization of the present work to the electrovac- 
uum case is important to establish the influence of the mag- 
netic dipole and quadrupole moments on the orbital motion of 
particles around compact objects (see e.g. Bakala et al. 2010; 
Sanabria-Gomez et al. 2010; Bakala et al. 2012). 

Interesting effects on the epicyclic frequencies due to the 
presence of the magnetic dipole have been already pointed 
out recently by Bakala et al. (2010) and Bakala et al. (2012). 
These effects were predicted after neglecting the contribu- 
tion of the electromagnetic field to the curvature, for y = 
see Bakala efaL (2010) and for j Bakala efaZ. (2012). 
Bakala et al. (2010) assumed the model of the star as a dipole 
magnetic field superposed on a Schwarzschild black hole. 
In the second work, they studied the case of a magnetized 
slowly rotating neutron stars; to build the model they super- 
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APPENDIX 
METRIC FUNCTIONS 

The functions A, B, C, H, G, K, and / used to express the metric functions (17) are given by 

^= X aijknrjn, ^Wirp C= ^ Cijnrj, ^ aijk, (Al) 

l<i<j<k<6 l<i<j<6 l<i<j<6 \<i<j<k<6 

H ^zA-{Pi+P2+l3i)B+ ^ hijkri rj + ^ (ai + aj) bij n rj, (A2) 

1<!<><A<6 l<i<y<6 

G^-(J3i+/32+/3i)A+zB+ gijrirj+ ^ (ai + aj + aOaijkri rj n, (A3) 

I</<y<6 \<i<j<k<(3 

6 

I ^if\+fi+ -B) + (/3i+P2+/33-z)C+ Y PUk^i n ru + ^ pi n + ^PU ' ('^'^ + o'v)'^-;]'"; fj' (A4) 

1 <i< j<k<6 != 1 1 <!< j<6 

with 

n = + (z- a,)2 , aijk = (-l)''^^'+'A,y,r,|,„„ , bij = {-irjAijHu^np , 
Cij = {-!)'''■' Aij[f{ai) T„,i„p - f(a,„) r„|p/ + f{a„) Tpii,„ - fia,,) T/i,,,,,] , 

hijk - (-i)'^'^'^^ijk(e*i 622\lmn + ^2 ^31|/mn + 6i2\lmn) , gij - (" 1 )'^^'^iv(a/ Tmlnp - a„, r„|p/ + ff,, Fpi/,,, - apTumn) , 

Pi = (-1)'D,[/(q'/) H,„\,y„ - f(a„) H„\p,i + f(a„) Hp\^,,,n - f(ap) H,^,,nn + fia.,) H,\,„„p] , 

Pi] - {-^y^^ Aij{e\ T23|/mnp + T3i|/,„„p + 63 Ti2\lmnp) , Pijk - (- 1)'^''^' A,jt(e j ^27,\lmn + e\ ^3\\lmn + ^3 ^ \2\lmn) , 

Aij = (a,- - aj) DjDj , A,;t = (c,- - ajXa/ - UkXaj - a^) D/Dj Dk , 

Di - — — — ■ , r/|„„, = Hiiai) Ai2|„m + H3{a,„) Ai2|,,/ + H3(a„) Am/,,, , 

(a; -PiKai - /32)iai - Pz) 
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and 



Hi(a„) Hjop) - H,{ap) H,„(a„) , Hi(a„) 



Z '^fl fk 
. . W* -Bi){a„- 



nL,os; - PI) nil (ys; - a) ti^'i- ^*)(^« - 



^lm\nps ~ ^im\np ^lm\ps ^lm\sn-> ^l\mnp ~ ^'ii.^l) ^\2\mnp ■> Hl\mnp ~ ^^l\mnp "t" ^^m\iipl "t" ^^n\plm ^" ^^p\lmn i 
^lm\nps - fia,,) Aii„\ps + fiOp) A;,„|,„, + /(a.,) Ai,„i„p , 

'^lm\nprs - f(a„) 6lm\prs - fi^p) Sl„,\rsji + fi^r) Slm\snp - /(ttj) Slm\npr , 

being a's the roots of the Sibgatullin equation SibgatuUin (1991); Manko and Sibgatullin (1993) 

e(z) + g(z) + 2/(z)/(z) = 0. (A5) 

KERR'S METRIC IN WEYL-PAPAPETROU QUASI-CYLINDRICAL COORDINATES 

In order to keep comparisons in the save place, we consider useful to display the Kerr solution in the Weyl-Papapetrou quasi- 
cylindrical coordinates. For this case. 



/ 



AA-BB 
(A-B)(A-B) 



„2r 



AA-BB lm[(A + B)H -(A + B)G] 



KKY[r„ 



AA - BB 



(Bl) 



where for our own convenience we do not present the definition of each term, but present the final combination of them, i.e., 



AA - = -8 (fl^ - m^)^ [p^ + z^) nr ^J-2z y/m--a^-a^ + + + ^/( Vm^ - + z) 



I ( Vm^ -a^ +z) +p2 

- 2fl- ^J -2z Vm^ _ ^2 _ q2 + ,„2 + p2 ^2 ^ ^ -^^2 + ^2^2 _ 

(A - B)(A - B) = -8 (m^ - a^f 



(B2) 



4 2 2 2 2 

m + m'p + m z 



— 



X (p^ + z^) |2m | -y/ -2z Vm- - - + m'^ + p^ + z~ + yjlz Vm^ - _ ^2 + ^2 p2 ^. ^2^ 
+ 2 ■sj—2z Vwi- - - + trP- + p^ + z~ sj 2z Vm^ - — -H -I- -I- -t- 3m^ j 

^ ■\j—2z Vm^ - fl^ — fl2 -t- -(- p2 -I- -H 2z Vm^ - - -1- -t- p^ -1- j 

H- OT^ -H p2 H- z^ -\/2z Vwz^ - - + m~ + p^ + z~ + 3m^ + p^ + jj , 
KkY] r„ = 16 (m^ - a^f (p^ + z^) ^^^(zW5^5)%p2 ^( Vm^ - ^2 -n z)% p^ 

! (m^ - fl^)"' (p^ 4- z^) ^-OT^ -\/ -2z Vm^ - - + + p^ + z^ 

— m ■\J—2z Vwi^ — fl2 - fl2 + ,,^2 + p2 ^ ^2 ^ Vot^ _ q2 ^. -)- p2 

H- -sj -2z Vwz^ — fl2 — + ,,^2 + p2 ^. ^2 _ ^ Vm^ — ■\J—2z Vm^ — _ ^2 ^. ,^2 + p2 ^. ^2 

- ( Vm^ -"fl2 H- 2^ -i-"p2 + q2 ^ ^ V„^2 _~a2 + 



(B3) 



(B4) 



lm[(A H- B)// - (A + B)G] = 16flm 



(B5) 



-H z Vot^ — ^ ^ Vm^ — + 7^ h- p^ -1- — m'' 4- mp^ + mz^ 



From here, it is clear how changing a -a will cause only a global change in the sign of the metric function u and therefore 
only a change in the g,^ metric component. 



